Evidence that solar flares drive global oscillations in the Sun 
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ABSTRACT 

Solar flares are large explosions on the Suns surface caused by a sudden release of magnetic 
energy. They are known to cause local short-lived oscillations travelling away from the explosion 
like water rings. Here we show that the energy in the solar acoustic spectrum is correlated with 
flares. This means that the flares drive global oscillations in the Sun in the same way that 
the entire Earth is set ringing for several weeks after a major earthquake like the December 
2004 Sumatra-Andaman Earthquake. The correlation between flares and energy in the acoustic 
spectrum of disk-integrated sunlight is stronger for high-frequency waves than for ordinary p- 
modes which are excited by the turbulence in the near surface convection zone immediately 
beneath the photosphere. 

Subject headings: Sun: flares — Sun: helioseismology 



1. Introduction 

The global 5-minute p-modes oscillations of the 
Sun have in the latest three decades provided us 
with valuable information about the physics in- 
side the Sun. The oscillations are believed to be 
excited by the turbulence in the near surface con- 
yection zone immed iately beneath the photosphere 
(jGough et al. 19961) . Complete reflection is ex- 



pected to take place only for oscillations with fre- 
quencies lower than the solar atmospheric acoustic 
cut-off (about 5.3 mHz). Oscillations with higher 
frequencies, which are seen in the Sun up to 9 mHz 
(sometimes known as high-frequency waves) are 
only partially reflected and have transmitted com- 
ponents, whi ch propagate upwards thr ough the 
atmosphere (|Balmforth fc Goughl ll99oh . While 
models of near surface convection can now ex- 
plain the excitation and damping of the ordinary 
p-modes no consensus has been reached about the 
details of what drives oscillations with frequen- 
cies larger than the acoustic cut-off frequency in 
the solar atmosphere. Different studies have pro- 
posed different models of the high-frequency waves 



Balmforth fc Goughl Il990t iKumar fc Lul Il991 
Jain fc Robertslll996l ). but none has explained all 
the features seen in the observations. 

We have used data from the SOHO (Solar and 



Heliospheric Observatory) (jFrohlich et al.l [1995) 
and the GOES (Geost ationary Ope rational En- 
vironmental Satellite) ( Garcia! Il994 ) satellites to 



show that there is a strong correlation between 
the energy at high frequency in the solar acous- 
tic spectrum and the appearance of solar flares. 
The correlation suggests that the driving of the 
high frequency waves is related to solar flares. 
Such a driving mechanism is well known from 
earthquakes where the entire Earth is set ring- 
ing for several weeks after a major earthquake 
like the December 2004 Sumatra- Andaman Earth- 
quake ([Park et al. I l2005h . In fact it was shown al- 
ready in 1972 that solar flares could excite free 
global oscillations in the Sun in the same way 
large eart hquakes exc ite free global oscillation in 
the Earth (|Woln1ll972[ ) . Since then a few attempts 
have been made in order to find a correlation be- 
tween the energy of the low degree p-modes and 
flar es but none of them have been very success- 
ful dGavrvusev fc Gavrvuseva l 19991 : Chaplin et al 



2004; I Ambastha fc Antiall2006f) . Instead flares are 



found to cause local short-lived high-degree oscil- 
lations in the active regions in which they are pro- 
duced ( Kosovichev fc Zharkova 19981 ). 

High-frequency waves were di scovered in high - 
degree observations of the Sun (jLibbrechtJ [l988) . 
but recently they have also been seen in disk in- 
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Fig. 1. — Smoothed power density spectrum of the 
Sun. The spectrum consists of granulation noise 
(below 1.8 mHz), acoustic power from completely 
trapped p-modes (between 1.8 and 5.3 mHz), and 
power from the high-frequency waves (with fre- 
quency higher than 5.3 mHz). The peak at 5.55 
mHz is an artifact. 



tegrat ed data such as the BISON radial velocity 
data (IChaplin et al.ll2003f) the GOLF radial ve- 
locity daia^JimiMz^t^^Ji^^aiid the VIRGO 
intensity data (j Jimenez et al.ll2005h . As the high- 
frequency waves have the highest S /N in intensity 
(j Jimenez et al.l [20051 ) we have chosen to focus on 
Virgo data (green channel) in this study. Fig. 1 
shows a power density spectrum of the solar data. 
Here the high-frequency waves appear clearly from 
5.3 mHz up to almost 8 mHz. 



2. Frequency/time diagrams 

We have calculated frequency/time diagrams of 
the integrated sunlight in order to evaluate the 
temporal behaviour of the high-frequency waves. 
This has been achieved by calculating power spec- 
tra of a number of running substrings of a long 
time series from the VIRGO instrument on SOHO. 
The pow er spectra we re calculated as least squares 
spectra ( Lombl 1976 ) in order to insure proper 
treatment of the gapes in the data. We tested 
different length of the substrings and found 7.5 
days to be the smallest length that we could use 
and still have enough S /N in the spectra to see the 
high-frequency waves clearly. In order to keep the 
rest of the analysis as simple as possible we calcu- 
lated the spectra with a frequency resolution of 1 



/iHz. The power spectra were converted to power 
density spectra by multiplying a given power spec- 
trum with the effective length of a given substring, 
which we calculated as the reciprocal of the area 
under the spectral window (in power) of the sub- 
string (jKieldsen et al.ll2005f) . 

Successive, overlapping, substrings were dis- 
placed by 18 hours with respect to each other. The 
procedure yielded 4620 power spectra. The fre- 
quency/time diagrams can now be constructed by 
stacking the power spectra vertically, producing a 
power diagram with frequency varying along the 
horizontal axis and time on the vertical axis. In 
order to enhance the visibility the diagrams were 
smoothed by a Gaussian point-spread-function 
having widths of 33 /iHz in frequency and 24 days 
in time. The widths were chosen to be as small as 
possible and still allow the structures of the high- 
frequency waves to appear as clear as possible. 

As the analyzed VIRGO data have a sampling 
of 60 sec and thus a Nyquist frequency of 8.3 mHz 
we have chosen to calculate the power spectra up 
to 8 mHz. There are some gaps in the data set, 
mainly around the SOHO vacation in the sum- 
mer of 1998, but when we smooth the diagram 
interpolation is made of these gaps by the Gaus- 
sian point-spread-function. As we generally do not 
think that this interpolation is reliable we have re- 
move the data in the gaps from Figs. 2 & 3. 

The frequency/time diagrams in Fig. 2 exhibit 
a regular pattern of vertical bright lines at low fre- 
quency. These indicate the power in the ordinary 
p-modes. But the pattern continues beyond 5.3 
mHz, demonstrating that underlying mode struc- 
ture is present also in the high-frequency waves. 
The frequency/time diagrams also show a pattern 
of horizontal lines. It is evident from the mid- 
dle panel that this pattern is correlated with so- 
lar flares. The number and strength of the flares 
has been represented by the X-ray flux measured 
in the soft channel (1- 8 A) of the X -ray sensors 
on the GOES satellites (|Garcialll99i) . The X-ray 
data have been smoothed with a Gaussian running 
mean of width 25 days in accord with the resolu- 
tion of the oscillation data. The width was taken 
to be 25 rather than 24 days as the Virgo data to 
some extend has been smoothed two times, first 
by calculating a power spectrum of a substring of 
length 7.5 days and second by smoothing with a 
Gaussian point-spread-function of width 24 days 
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Fig. 2. — Smoothed frequency /time diagram of the Sun from disc-integrated data from the VIRGO instru- 
ment on SOHO is shown in the left panel. The clear vertical lines in the left of the image are the p-modc 
oscillations, and it is evident that these lines continue far beyond 5.3 mHz. The vertical line at 5.55 mHz 
is an instrumental artifact. The bright horizontal lines are strongly correlated with solar flares, represented 
by the X-ray flux in the middle panel. These lines are most prominent in the high-frequency region of the 
spectrum. The right panel shows the amplitude spectra of the individual power density spectra in the left 
panel, with a bright line at the frequency separation of the high-frequency waves of 0.014 /iHz -1 (70 /iHz). 
The gray regions mark times with no data from SOHO. The color scale is logarithmic. 



in the vertical direction. 

3. Correlation between solar flares and 
power density at different frequencies 

In order to analyze how flares correlate with 
oscillations at different frequencies we have com- 
pared the variation of the total power within 20 
fiHz of three different modes (at 3.22, 5.24 and 
7.22 mHz). It is seen in Fig. 3 that the correlation 
is strongest for the high-frequency oscillations at 
7.22 mHz, at 5.24 the correlation is much weaker 
and at 3.22 the dominant feat ure is the well known 
slow variation of the dynamo (jChaplin et al. 2003. 
and references herein). All the three modes have 
higher power density around 2006 than around 



1996, though the power is expected to be the same 
as both times are close to solar minimum. This 
trend is also seen in the X-ray flux (the mean 
power density is higher between 2004 and 2006 
than between 1996 and 1998), but as the trend 
is monochromic increasing it could be an instru- 
mental effect. Though we do not know of such an 
effect in the instrument. 

In Fig. 4 we show the general correlation be- 
tween the power density and X-ray flux indepen- 
dent of the individual modes. In order not to 
have the correlation dominated by the long pe- 
riod trends we have subtracted a Gaussian run- 
ning mean with a width of 2 years from both the 
power density and the X-ray flux. Fig. 4 shows 
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Fig. 3. — A zoom in on Fig. 2 where the regularity 
of the high-frequency waves appears more clearly. 

that the temporal behavior of the power density 
spectrum is more correlated with the temporal be- 
havior of X-ray flares at high-frequency than at 
low frequency. 

We do not know the exact excitation mecha- 
nism for p-modes being excited by flares, however 
two possible effects may be considered in order to 
explain and understand the behavior of the high- 
frequency p-mode power. First of all the observed 
power of p-modes at high-frequency are signifi- 
cantly lower than the power pe r mode observed at 
the pe ak (3 mHz). According to lStein fc Nordlund 
(|200ll ) one find the p-mode power to decrease in- 
versely to the fourth power of the frequency (above 
3 mHz) simply as a consequence of the granula- 
tion power decreasing as i/~ 4 . Therefore a sin- 
gle event (like a flare) that excites a mode will 
have a larger relative effect at high frequencies 
because the other excitation sources are much 
smaller. Related to this one should also note 
that the high-frequency background (between p- 
modes) show an increased correlation with flare 
activity. One may therefore argue that the back- 
ground noise from granulation correlate with flare 
activity (as clearly seen in Figs. 2 and 5) and the 
increased background transfer more power into the 
high-frequency modes via the stochastic excita- 
tion, which is exactly what is observed. At fre- 
quencies around 3 mHz one do not find the back- 
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Fig. 4. — The X-ray flux from the Sun (lower 
panel), total power of the mode at 7.22 mHz (sec- 
ond panel), at 5.24 mHz (third panel) and at 3.22 
mHz (top panel). It is evident that the correla- 
tion between the total power of the modes and 
the X-ray flux increases with increasing frequency. 
The slow change of the total power of the mode 
at 3.22 mHz is a result for the solar dynamo 
(jChaplin et al.ll2003l . and references herein). 



ground to correlate as strongly (see Fig. 5) as at 
higher frequencies (5-7 mHz). 

This may explains why attempts to see a cor- 
relation between the energy of the ordinary p- 
mo des and flares have not been very success- 
ful (iGavrvusev fc Gavrvusevalll999l ; Chaplin et al 



2004 lAmbastha fc Antial 120061 ) . Fig. 4 shows a 
steady clime of correlation over the acoustic cut- 
off frequency, not a sudden jump. This is in agree- 
ment with what we expect as the reflectivity of the 
solar atmosphere is expected to decrease continu- 
ously over the acoustic cut-off frequency. 
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Fig. 5. — Correlation between the power density 
and the X-ray flux from the Sun as a function of 
frequency. 



4. Correlation between solar flares and 
global oscillations at high frequency 

The acoustic spectrum contains two parts at 
high frequency: a background and the acoustic 
waves. In order to separate these two parts we 
ha ve preformed the sa me analysis as preformed 
by I Jimenez et al. ( 2005 ) and calculated amplitude 
spectra of the individual power density spectra in 
the frequency/time diagram. The exponential de- 
cay of the power density spectra were removed be- 



fore thes e amplitude spectra were calculated as 
done by iJimenez et alj (2005). The amplitude 



spectra as function of time are shown in the right 
panel of Fig. 2. It is seen that the frequency sep- 
aration of the high-frequency waves remains con- 
stant at 70 /iHz throughout the observation inde- 
pendent of solar flares. 

5. Discussion 

The correlation we see between the power den- 
sity of high-frequency wave and the flares is similar 
to what one would expect from a comet impact. 
This case was analyz ed for the Sh oemaker-Levy 
impart on Jupiter by I Gough ( 1994 ). In that pa- 
per it was also noted that such impacts on the 
Sun would have the largest response in the high- 
frequency domain. The discovery presented here 
will improve our understanding of the flares. An 
advance here will be provided by the asteroseismic 
data that we shall receiv e from the Kepler satel- 
lite ( Borucki et al.l 120031 ). High-frequency waves 



have been detected in asteroseismic data of the 
solar-like sta rs beta Hydri and alpha Cen A & B 
(Karoff 2007), and with the Kepler data we should 
be able to observe the time variation of the ampli- 
tudes of such waves in a large number of solar-like 
stars, perhaps suggesting the presence of flares in 
those stars. The last outcome could aid our un- 
derstanding of stellar dynamos. 
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